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Abstract. This paper is the first of a series dedicated to the X-ray properties of the young open cluster NGC 6231. Our data set 
relies on an XMM-Newton campaign of a nominal duration of 180 ks and reveals that NGC 623 1 is very rich in the X-ray domain 
too. Indeed, 610 X-ray sources are detected in the present field of view, centered on the cluster core. The limiting sensitivity 
of our survey is approximately 6 x 10~ 15 ergcirT 2 s _I but clearly depends on the location in the field of view and on the source 
spectrum. Using different existing catalogues, over 85% of the X-ray sources could be associated with at least one optical and/or 
infrared counterpart within a limited cross-correlation radius of 3" at maximum. The surface density distribution of the X-ray 
sources presents a slight N-S elongation. Once corrected for the spatial sensitivity variation of the EPIC instruments, the radial 
profile of the source surface density is well described by a King profile with a central density of about 8 sources per arcmin 2 and 
a core radius close to 3.1 arcmin. The distribution of the X-ray sources seems closely related to the optical source distribution. 
The expected number of foreground and background sources should represent about 9% of the detected sources, thus strongly 
suggesting that most of the observed X-ray emitters are physically belonging to NGC 623 1 . Finally, beside a few bright but soft 
objects - corresponding to the early-type stars of the cluster - most of the sources are relatively faint (~ 5 x 10~ 15 erg cirT 2 s ) 
with an energy distribution peaked around 1.0-2.0keV. 

Key words. Open clusters and associations: individual: NGC 6231 - X-rays: individuals: NGC 6231 - X-rays: stars - Stars: 
early-type - Catalogues 



1. Introduction 

Detailed studies of young clusters are powerful tools to probe 
crucial astrophysical issues. Because they a priori contain both 
early-type stars and pre-main sequence (PMS) stars, young 
clusters are privileged laboratories to test star formation and 
evolution theories. They indeed provide a homogeneous sam- 
ple of stars in terms of distance, reddening, environment, chem- 
ical composition and age. With the currently available X- 
ray observatories, unprecedented investigations of young open 
clusters in the X-ray domain have been performed in the past 
few years. The increased sensitivity, spectral power and reso- 
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lution of the XMM-Newton and Chandra observatories, com- 
pared to X-ray satellites of the previous generations, give now 
a much more complete view of the X-ray properties of the star 
populations in clusters. 

For example, a 76 ks Chandra observation of the embedded 
young cluster NGC 2024 (d ~ 410 pc; age = 0.3 - a few Myr) 
revealed 283 X-ray sources displaying heavily absorbed hard 
spect ra with a mean temperature kT ~ 3 ke V ( Ski nner et alJ 
l2003h . A significant fraction (25%) of the X-ray sources shows 
a wide range of variability within the exposure duration. In 
addition, Chandra detected at least 96% of the known classi- 
cal T Tauri st ars in NGC 2024 . Results for other clusters are 
very similar. iRauw et all J2003I) performed a 20 ks observation 
of NGC 6383 (d ~ 1.4kpc; age = 1.7 - 5 Myr) and found 77 
sources, mostly centered on the cluster location. An important 
fraction of these sources are probable PMS obje cts. Using both 
XMM-Newton and Chandra facilities, Preibis ch & ZinneckeJ 
ll2004 and references therein) studied the very young stel- 
lar cluster IC348 (d ~ 310pc; age ~ 2 Myr) and found 286 
X-ray sources among which over 50 classical T Tauri stars. 
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Comparison of Chandra- and XMM-Newton- based spectral 
properties suggested that the X-ray characteristics of T Tauri 
stars remain mostly constant over periods of years. NGC 6530 
(d ~ 1.8 kpc; age ~ 1.5-2.0Myr) is a very rich open clus- 
ter containing several massive O-type stars as well as a large 
population of B-t ype stars. XMM-Newton observations by 
iRauw et all J2002) revealed 119 sour ces, of which a large frac- 
tion are PMS candidates. Similarly to lSkinner et alJ results, the 
X-ray spectra of the PMS candidates are characterized by tem- 
peratures of a few keV. Using a 60 ks Chandra observation 
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centered on the same cluster, Damiani et al. (2004) revealed 

884 X-ray sources, among which 90 to 95% are PMS stars. 

From this review of the recent literature, there is an obvi- 
ous body of observations showing that, besides the expected 
X-ray emission from massive stars, a large population of X- 
ray emitting low-mass PMS stars is to be found while observ- 
ing young clusters. The present study of the very rich cluster 
NGC 623 1 lies in this framework. It aims at a better compre- 
hension of both early-type stars and young open clusters by ex- 
tending the previous sample of investigations to clusters with 
a large O-type star population. A severe limitation of several 
of the above cited works is the lack of detailed studies on the 
concerned cluster at other wavelengths. Indeed, with Chandra 
and XMM-Newton, the X-ray observations are so deep that a 
deep optical photometry of the field is required. Such a data 
set is indeed essential to, for example, more clearly identify the 
evolutionary status of the different sub-populations of the clus- 
ter. Fortunately, as shown by the literature review of Sect. [3] 
the stars in NGC 623 1 have been thoroughly studied. Together 
with the depth of the present X-ray campaign, this is one of the 
strengths of the current work. Finally, the present work dis- 
tinguishes itself from the previous investigations because of 
the particular planning of the X-ray observations. Indeed our 
XMM-Newton campaign towards NGC 623 1 was actually split 
into six successive pointings, spread over a period of five days. 
This will allow us to probe the variability of the X-ray emission 
of the detected sources on different time-scales. 

A detailed analysis of the central target of the field, the col- 
liding wind binary HD 152248, has been presented recently in 
a dedicated paper (San a et al.l2004l) . The source will therefore 
not be further discussed in more details in the present paper. 
Preliminary results from this campaign, mainly rel ated to the 
early-type X-ray emitters, were also presented in ISana et all 
(2005b, 2006b c). In the present paper, we focus on the X- 
ray catalogue and we discuss some general properties of the 
detected sources. Other aspects of the X-ray properties of 
NGC 6231, such as the early-type and the pre-main sequence 
population characteristics, will be addressed in subsequent pa- 
pers of this series. 

This first paper is organized as follows. After a review of 
the abundant literature on NGC 623 1 and on the Sco OB 1 as- 
sociation, Sect. [3] describes the campaign and the subsequent 
data reduction processes. In Sect. @J we address the detection 
and identification of the sources in the XMM-Newton field of 
view, and we present the resulting X-ray catalogues. Finally, we 
probe the main properties of the cluster X-ray emitters (Sect. 
EJ. Sect.[6]summarizes the results of the present work. 



2. NGC 6231 and the Sco OB 1 association: a 
literature review 

2.1. The Sco OB 1 association 

Located in the Sagittarius-Carina spiral arm of our galaxy 
(g(2000) = 16 h 53 ?6, <5(2000) = -41°57'; I = 343!3, b = 1°2, 
|PerrvetalJfl99l . the Sco OB 1 association is an extremely 
rich and interesting region of the sky. 2° long by 1° wide, 
it extends from the gaseous nebula IC4628 on its northern 
end to the young open cluster NGC 623 1 towards its southern 
end. Its majo r axis is appr oximately parallel to the Galactic 
plane l lMorgan et alJll95 3a). A sparser group, Tr24, is to be 
found near IC 4628 while two other clusters, NGC 6242 and 
NGC 6268, lie slightly north of the association. Finally the H n 
region IC4878, centered on NGC 6231, extends by about 4° 
by 5° in the form of an elliptical ring and is probably triggered 
by the cluster. The emission nebula is faint within the ring but 
is very bright where the ring is crossed by the northern end of 
the association dBok et al.ll966h . 

The interest in Sco OB 1 mainly originates from its ex- 
tende d early-type star content (Shobbrook 1983: lRaboud et all 
Il997l) . Beyond the numerous O- and B-type stars, the associa- 
tion also shelters two of the rare Wolf-Rayet (WR) stars, two 
Of stars di splaying P Cygni profiles as well as several B Cephe i 
variables dBalona & Engelbrechtlll985t lArentoft et al1l200lh . 
Among the peculiar objects found within Sco OB 1 is the bright 
star Sco. With an absolute magnitude around My=-8.3, 
( 1 Sco is one of the brightest stars of the Milky Way. Many of 
the 'normal' early-type stars further present signs of variability 
and have a good chance to be binary systems ("e.g. iRaboud 
1996[ lArentoft et alJ200lT) . 



2.2. The NGC 6231 cluster 

Located near the southern end of the Sco OB 1 association, 
the young open cluster NGC 6231 (ff(2000) = 16 h 54 m 09 s , 
£(2000) = -41°49'36") contains an important number of 
bright early-type stars in its centre. Often c onsidered as the 
nucleus of the association dBoketalJll966l) . its relationship 
to Sco OB 1 h as been subject to differen t interpretations with 
time. Though iHeske & Wendke^ dl984l) presented evidence 
that Tr24 and Sco OB 1 form a single aggregate, these authors 
proposed that NG C 6231 is actually a foreground cluster. 
IHeske & WendkeJ also found a sub-cluster of PMS s tars in the 
vicinity of Tr 24. Based on an extensive set of data. lPerrv et alJ 
dl990l 1991) re-addressed these issues and carefully studied 
the interrelation between the three aggregates. They estab- 
lished that Sco OB 1, NGC 6231 and Tr24 are located at the 
same distance and have the same age, thus demonstrating 
that NGC 6231 is not a foreground object but is clearly 
embedded in the Sco OB 1 association. NGC 623 1 therefore 
retains its status as the nucleus of the association. IPerrv et all 
could ho wever not confirm the three stellar sub-aggregates 
found by Seggewis 1 fl968a) in Tr24 and, as suggested by 
IHeske & Wendkerl dl984h . "thev casted further doubts on 
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the physical reality of the Tr24 aggregate itself. IPerrv et alJ 
finally confirmed the existence of a PMS sub-cluster near Tr 24. 

The properties of NGC 6231 and of its stellar con- 
tent have been thoroughly investigated during the past 
century. Three main streams of investigation were indeed 
designed, namely photometry, spectral classification and 
radial velocity measurements. The photometric approach is 
however predominant and was extensively performed using 
different photometric systems. The bulk of the available 
literature on the cluster re lies on photograp h ic, ph otoelec- 
tric or CCD campaigns: Brownlee & Cox ll 19531 PV), 
Houckl dl956h. iWair ven & Walraven ( 1960, Walraven), 



I dl956l). iWalraven & Walravenl d!96ol 
Feast et al J (1 1 96 ll UB VI. iBreckinridee & KrorJ d 1 9631 PV), 
; all ill 9661 UBV Hffl. lFeinstein & Ferrerldl968l UBV). 



Bok et al. ( 1966, UBV H/3), Feinstein & Ferrer 
Seggewissl dl968bl UBV). ISchild et al] ill 9691 UBV), 
Crawford et alJ dl97ll nvhv HB). Ic.arrison & Schildl dl979| 
UBVl lShobbrookldl983l uvbyH/3) . lHeske & Wendkerld 19841 
UBV) Ivan Gende ren et alJ dl984l Walra ven), IPerrvetal] 
199ll uvbvl. iMevnet et all dl993l U BV). IBalona & Lanevj 
u vby HB). IRaboud et alJ d 19971 Gen eva). ISung et alJ 
UBVrRIV Ha). iBaume etaTT dl999L UBVI C ). The 
more recent works (from ~1990's) offer a much more complete 
view of the cluster both in terms of their angular extent and 
of the magnitude limit reached. An extensive still careful 
comparison of most (if not all) the works pub lished prior to 
1990 has been performed bv lPerrv etal.ldl99ll) . 



1995 



1998 



Spectral classificat ion of the cluster ob j ects has mainly 
been carrie d out by Morgan et al. 1 dl953bl). iHouckl fl 956). 
iFeast etall dl96ll). ISchild let alJ (M969I). iGarrison & Schildl 
dl979l).lLevato & Malarodaldl980l) Iconti & Alschulerldl97ll) 
IWalbornl Jl972l). iMathvsl Jl988l Il989h . Idarcia & Mermi lliod 

Ra dial v e locity campaigns 



i vvaiuyiu LYiauiva ill 

(2001) and ISanal |2005 



iLevato & Morrelll i 


1983), Levatoetal. 


dl988l). 


Perry et alJ 


(1990), Penny etal. 


(1994), Stickland & Lloyd (2001, IUE 


data). 


Garcia & Mermilliod (2001). LSanaetal.1 d2002k and 



Sana (2005). Several authors also paid a special attention 
to particular objects, mainly binaries of which they per- 
formed a more detailed study. These o bjects are WR79 



dLuhrslll997l), HP 152 2 18 dStickland et al.|ll997l). HP 1 52248 




( Stickland et al. 
ISanaetal.1 12004 



19961 IPennv et al.l 119991 ISanaetafl 12001 



), 



CPP-41°7742 
and HP 152219 (|Sana^LaJ2006a). 



boo it 

2005a) 



Aside from these three main streams, several authors 
addressed specific aspects of the cluster that provide a use- 
ful complementary view. Among other topics, photometric 
variabil i ty of a few dozens of objec t s was investigated by 
iBalonal dl983h . IBalona & Engelbrechfl 1 1985 1. IBalonal ill 992) 
and more recently by lArentoft et al.1 pOOl). These studies 
allowed to detect several ySCephei, a couple of 6 Scuti and 
a few other variable stars, including a coupl e of eclipsing 
binaries. The binary fraction was estimated by[Raboud ( 1996) 
and lGarcia & M ermilliodl feoOll) . iRaboudl derived a minimum 
binary frequency of 52% in his sample o f 53 B-type stars with 
a spectral type between Bl and B9 while iGarcia & Mermilliodl 



(2001) obtained an extremely high frequency of 82% for stars 
earlier than B1.5V an d, in particular, of 79% for th e O-type 
stars of the cluster. IRaboud & Mermilliod (1998) showed 
evidence of mass segregation in NGC 6231, most probably 
related to the formation processes rather than to the dynamical 
evolut ion o f the c l uster. Proper motions were studied bv lBraesI 
( 1967J) and lLavall dl972) while most of the O-type stars of the 
clusters were includ ed in the large ICCP Speckle campaign of 
iMason et al.ldl998l) . 



The distance modulus (DM) of th e cluster reported in 
the e arlier literatur e ranges from 10.7 ( Merm illiodl 1 1 98 1() to 
11.9 dHoucklll956l 2300pc - cited b y iBok et alJll966l) . In a 
more recent work, Per ry et al.l 11991) obtained DM = 11.50 
and 11.55 for ScoOB 1 and NGC 623 1 respectively , with an 
uncertainty of about 0.32. IBalona & Lanevl dl995l) derive d 
DM = 11 08 + 0.05 for NG C 6231: IRaboud etall d 19971) 
1 1 .2 ± 0.4: ISun g et alJ dl998l) . 1 1 .0 + 0.07 and IBaume etall 
(1999) 11.5 ± 0.25. The weighted mean of these five latter 
values gives DM = 11.07 ± 0.04, corresponding to an actua l 
distance of 1637 + 30 pc. The same authors (but ISung et al J) 
respectively derived ages of 7.9 ± 0.9 Myr, 5 + 1 Myr, 
3.8 + 0.6 Myr and 3 to 5 Myr. On the basis of the R-Ha index, 
ISun g et al.ldl998l) found 12 PMS stars plus 7 PMS candidates. 

A controversial question is the probable differential redden- 
ing across the cluster. Such a differential reddening was first 
suggested by Breckinridge & Kron (1963), outlining that the 
southern part of the cluster suffers a heavier reddening. O ther 
autho rs dShobbrooklll983l IPerrvet aljfl99ll IBalona & Lanevl 
1995) rather pr oposed a unif o rm red deni ng across th e field. 
More recently, Raboud et ail d 19971) a nd ISung et alJ dl998l) 
resul ts strongl y support the first idea of iBreckinridge & KrorJ 
(1963J), and ISung et al.l presented a map of the reddening 
distribution in NGC 623 1 . There seems to be an agreement in 
the early literature that most of the reddening occurs between 
a distance of 100 and 1300 pc. Based on FUSE observations, 
Marggraf et ID §004) recently confirmed angular variations 
in the column density towards the core of the cluster. They 
reported that the absorption towards NGC 6231 occurs in 
several foreground clouds. The main absorption component 
lies in the Lupus cloud region at a distance of 150 pc, while 
the second one is probably in the vicinity of the Sco OB 1 shell 
surrounding NGC 623 1 . Finally, ICrawfordl d200ll) probed the 
structure of the interstellar Nai and Ki towa rds the cluster 
and reached conc lusions similar to those of Marggraf et alJ 
(2004). ICrawfordl also outlined that no clues of active shocks 
in the shell components could be foun d. Pol a rimetr ic ob- 
servations were performed by Feinstein et alJ (120031) who 
found evidence for a past supernova explosion in the cluster. 
These authors however suggested that their observations could 
also be explained by a bubble triggered by winds from hot stars. 

Finally, we note that the Hipparcos parallaxes derived for 
NGC 6231 were known to be problemat ic with a negative 
mean value of -0.8 + 0.4 mas ( Arenou & Luri 1999). These re- 
sults were recently revised by Makarov (2003a) who obtained 
1 .7 ± 0.4 mas, corresponding to a distance modulus of 8.9 ± 0.5, 
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Table 1. Journal of the XMM-Newton observations of NGC 6231. Columns 2 and 3 give the spacecraft revolution number and 
the observation ID. The Julian Date (JD) at mid-exposure is reported in Col. 4. Cols. 5 to 7 (resp. Cols. 8 to 10) list the performed 
(resp. effective - i.e. corrected for the high background periods) exposure times for the EPIC MOS1, EPIC MOS2 and EPIC pn 
instruments. The last column provides the position angle (PA) associated to the revolution. The total campaign duration is given 
in the last line of the table. 



Obs.# 


Sat. 


Exposure 


JD 


Performed duration (ks) 


Effective duration (ks) 


PA 




Rev. 


ID 


JD-2450000 


MOS1 


MOS2 


pn 


MOS1 


MOS2 


pn 


DDD:AM:AS.s 


1 


319 


0109490101 


2158.214 


33.3 


33.3 


30.7 


33.1 


33.2 


30.6 


274:57:11.5 


2 


319 


0109490201 


2158.931 


22.1 


22.1 


20.2 


19.8 


19.8 


16.5 


274:57:11.5 


3 


320 


0109490301 


2159.796 


34.4 


34.4 


31.8 


33.7 


33.9 


30.1 


275:35:26.6 


4 


320 


0109490401 


2160.925 


31.4 


31.4 


29.1 


26.0 


24.3 


11.7 


275:35:26.6 


5 


321 


0109490501 


2161.774 


31.1 


31.1 


28.5 


30.9 


31.0 


28.4 


276:13:34.9 


6 


321 


0109490601 


2162.726 


32.9 


32.9 


30.3 


32.9 


32.8 


30.3 


276:13:34.9 


Total duration (ks) 




185.2 


185.2 


170.6 


176.5 


175.0 


147.5 





however still far from the mean value obtained from the photo- 
metric studies. 

Turning to the X-ray domain, NGC 6231 was observed 
by the ROentgen SATellite (ROSAT). Thirty-five objects 
were detect ed, mainly associated wi th the early-type stars of 
the cluster. ICorcoran ill 9961 Il999h presented some results 
of this campaign as well as the X-ray light curve of three 
objects, namely HD 152218, HD 152248 and HD 152249. 
Only HD 152248 displayed clear variations of its flux. 
Finally a few objects were also o bserved at radio wavelength 
JSetia Gunawan et al jEooi F2003 ) but only half of them were 
detected. 



3. Observations and Data Reduction 

3.1. The XMM-Newton campaign 

The XMM-Newton campaig n t oward s NGC 623 1 has already 
been described in San a et alJ (120041 . For the sake of com- 
pleteness, we again give here a brief description of the X- 
ray observations. In September 2001, during sat ellite revolu- 
tions 319 to 321, the XMM-Newton observatory Jjansen et alJ 
12001 1) performed six successive exposures of an approximate 
duration of 30 ks. The field of view (FOV) was centered on 
the 07.5III+07III colliding wind bi nary HD 152248 ( a 2 <m = 
16 h 54 m 10!06, 5 2 ooo = -41 °49 , 30'/1 : ISana et aljEooTh . in the 
core of the cluster. Position angles (PAs) were very simi- 
lar through the six exposures, ranging approxi matively from 
274?95 to 276!23 All t hree EPIC instruments JStruder et alJ 
l200lUTurneretal.l 2001 ) were operated in the Full Frame mode 
together with the Th ick Filter to reject UV /optical light. The 
RGS spectrographs Jden Herder et alJboOll) were run in the 
Standard Spectroscopic mode. Due to the brig htness of the 
cluste r objects in the FOV, the Optical Monitor ( Mason et alJ 
l200ll) was switched off throughout the campaign. Table ^pro- 
vides the journal of the X-ray observations. 

3.2. Data Reduction 

The EPIC Observation Data Files (ODFs) were processed us- 
ing the XMM-Science Analysis System (SAS) v 5.4.1 imple- 



mented on our computers in Liege. We applied the emproc and 
epproc pipeline chains respectively to the MOS and pn raw data 
to generate proper event list files. No indication of pile-up was 
found in the data. We then only considered events with patterns 
0-12 (resp. 0-4) for MOS (resp. pn) instruments and we ap- 
plied the filtering criterion XMMEA_EM (resp. FLAG = 0) as 
recommended by the Science Operation Centre (SOC) techni- 
cal note XMM-PS-TN-43 v3.0. For each pointing, we rejected 
periods affected by soft proton flares. For this purpose, we built 
light curves at energies above 10 keV 1 and discarded high back- 
ground observing periods on the basis of an empirically derived 
threshold (adopted as equal to 0.2 and 1.0 cnts 1 for the MOS 
and pn instruments respectively). The so-defined GTIs (Good 
Time Intervals) were used to produce adequate X-ray event lists 
for each pointing from which we extracted images using x- and 
y-image bin sizes of 50 virtual pixels 2 . 

We finally combined the event lists obtained for all six 
pointings to increase the statistics of faint sources. For this pur- 
pose, we used the SAS task merge. For each EPIC instrument, 
we included the event lists resulting from different pointings 
one by one. We also built merged event lists that combine the 
twelve MOS or the eighteen EPIC event lists. The Attitude 
Files generated by the pipeline were merged using the same 
approach and we adopted, for handling the merged event lists, 
the Calibration Index File (CIF) and the ODF corresponding 
to the first pointing (Obs. 1 in TableQ. 



4. X-ray source detection and identification 

In this section, we focus on the detection and identification 
of the X-ray sources in the XMM-Newton FOV. For this pur- 
pose, we only used the merged event lists and images, ac- 
counting in this way for the six pointings at once. The total 
effective exposure times towards the cluster are, respectively 

1 Expressed in Pulse Invariant (PI) channel numbers and consider- 
ing that 1 PI channel approximately corresponds to 1 eV, the adopted 
criterion is actually PI > 10000. 

2 Though the physical pixels of the EPIC MOS and pn detectors 
have an extent on the sky of respectively 1'.' 1 and 4'.' 1 , the virtual pixels 
of the three instruments correspond to an extent 0705. The obtained 
images have thus a pixel size of 275. 
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Fig. 1. Combined EPIC MOS three-color X-ray image of the young open cluster NGC 623 1 . The field is roughly 30' in diameter 
and is centered on HD 152248. North is up and East is to the left. The different colors correspond to different energy ranges: red: 
0.5-1.0keV; green : 1.0-2.5keV; blue : 2.5-10.0keV. 



for the MOS1, MOS2 and pn instruments, of 176.5, 175.0 
and 147.5 ks. Together with the high sensitivity of the XMM- 
Newton observatory, the combination of the six pointings and 
of the three instruments provides one of the deepest X-ray 
views of a young open cluster. Fig.[^shows a three-colour im- 
age of NGC 623 1 and reveals a densely populated field with 
hundreds of point-like X-ray sources. This section therefore 
aims at providing a uniform catalogue of these sources. It is 
organised as follows. First we present the source detection pro- 
cedure as well as a brief description of the obtained catalogues. 



As a next step, we focus on the identification of the X-ray 
sources and, finally, we investigate the detection limit of the 
present data set. 

4.1. Source Detection 

We based our source detection on the SAS detection chain ede- 
tect_chain. For this purpose, we selected three energy ranges, 
a soft (S x ) band (0.5-1.0keV), a medium (M x ) band (1.0- 
2.5 keV) and a hard (H x ) band (2.5-10.0keV), and we built the 
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corresponding input images for the different instruments. The 
edetectjzhain task is formed by the succession of the S AS tasks 
eexpmap, emask, eboxdetect run in local mode, esplinemap, 
again eboxdetect run in map mode and finally emldetect: 

- eexpmap calculates the exposure maps corresponding to the 
input images; 

- these exposure maps are used by emask to build masks 
which select the relevant image areas where the detection 
should take place; 

- eboxdetect, run in local mode, uses a 5x5 pixel box and 
a surrounding background area to search for significant 
sources simultaneously in all input images; 

- esplinemap uses the resulting source list to remove sources 
from the input images and creates smooth background 
maps by fitting a 2-D spline to the source-subtracted im- 
ages; 

- run in map mode, eboxdetect uses a 5x5 pixel box and the 
values from the background maps to search for significant 
sources simultaneously in all input images; 

- emldetect finally uses the preliminary source list from 
eboxdetect and determines the source parameters (e.g. co- 
ordinates, count rates, hardness ratios, etc.) by means of 
simultaneous maximum likelihood psf (point spread func- 
tion) fitting to the source count distribution in all energy 
bands of each EPIC instrument. It also provides an equiva- 
lent logarithmic likelihood L2 (Ea. lA.U commonly used as 
an indication of the reality of the corresponding source. 

From our experience, the eboxdetect task run in map mode 
tended to eliminate some apparently real sources from the 
intermediate source list. We therefore preferred to use the 
preliminary source list obtained by eboxdetect in local mode as 
an input list for the psf fitting step performed by the emldetect 
task. This approach does not bias the result since, if the 
source is real, the psf fitting will provide a large logarithmic 
likelihood while, if instead the source is fake, the logarithmic 
likelihood will be low and the source will be rejected. Though 
more expensive in computation time, this approach results in a 
more complete source list. As it was known that the equivalent 
logarithmic likelihood values (L2) computed by the emldetect 
task in the SAS v 5.4.1 (and earlier versions) were erroneous, 
we implemented a patch to recover the correct L2 values. We 
give a brief description of it in Appendix lAl The problem is 
now fixed from SAS version 6.0 on. We checked our corrected 
logarithmic likelihood values against SAS v 6.0 results and 
found them to be in close agreement. 

We first performed single psf fit detection but, due to the 
crowdedness of the field, we also allowed for simultaneous 
fitting of up to four sources. In doing so, we adopted a value 
of 0.68 for the two parameters scut and ecut. This choice was 
led by the need to account for as large an energy fraction 
of the psf as possible while keeping the computation time 
down to reasonable limits. Due to the densely populated field, 
the wings of the source psf are often largely contaminated 
by emission from neighbouring sources. The adopted values 
therefore appeared as a reasonable compromise. On the axis, 
this corresponds to a physical radius of about 15". Only a few 



Table 2. Adopted detection thresholds for the equivalent log- 
arithmic likelihood L2 corresponding to the different EPIC in- 
struments (left column) or to any combination of them (right 
column). AppendixIBlprovides more details on how these val- 
ues were computed. 



Instr. 


Li 


Instr. Comb. 


L 2 


MOS1 


11 


MOS1 +MOS2 


21 


MOS2 


11 


MOS1 +pn 


35 


pn 


25 


MOS2 + pn 


35 






MOS1 + MOS2 + pn 


45 



tens of sources actually required multi-psf fitting, with three 
psf being simultaneously adjusted at the maximum. Finally, 
we re-ran the emldetect task allowing for extended sources to 
be fitted. A careful comparison of the resulting lists shows that 
only a few sources increase significantly their detection likeli- 
hood while allowing for extended source fitting. An inspection 
of the X-ray images and of the optical catalogues reveals 
that these sources most probably correspond to unresolved 
point-like sources rather than to physically extended sources. 

The described detection procedure was applied for each 
EPIC instrument as well as for any combination of them. The 
resulting source lists were generally consistent. The main dif- 
ference comes from the presence of different gaps in the differ- 
ent data sets. We built our final source list adopting the follow- 
ing criteria. 

(i) We selected the deepest combination of EPIC instru- 
ments, requiring that the detected source is distant by at least 
13" from any gap, bad column or detector edge. 

(ii) By a visual inspection of each source in all images 
and subsequent combinations, we adopted an empirical equiva- 
lent logarithmic likelihood (L2) threshold as the detection limit. 
This led us to consider the way to perform a consistent choice 
for the threshold values while dealing with different instrument 
combinations, and hence with different numbers of input im- 
ages. As a general comment, it is obvious that adopting a con- 
stant logarithmic likelihood threshold while dealing with differ- 
ent combinations of the EPIC instruments does not allow us to 
keep a constant threshold in terms of the signal level. Indeed, all 
other things being equal, the signal-to-noise ratio is increased 
while combining several detectors, allowing us in principle to 
detect fainter sources. However, in such a crowded field as the 
current one, we note that no important gain is achieved in terms 
of source detection. In other words, the very large majority of 
the detected sources are already seen in a single instrument, 
though of course combining the different instruments yields a 
much better estimation of their X-ray parameters. 

As a consequence, we have decided to adopt a logarith- 
mic likelihood threshold in one instrument and to look for the 
equivalent thresholds in any EPIC combination. This issue is 
presented into more details in Appendix 151 Table |2] gives the 
logarithmic likelihood thresholds finally used for the source 
detection. These values provide thresholds in various combina- 
tions that are consistent with the logarithmic likelihood-based 
detection threshold adopted in a single MOS instrument. We 



Table 3. Sample of the NGC 6231 X-ray source catalogue. Each part of the table refers to a particular EPIC instrument. First column gives the source number. Second column 
provides the source name (based on position in J2000.0) following the format specified for XMM-Newton sources. Note that the coordinates are truncated, not rounded. The 
estimated error on the 2D position (in ") is listed in the third column (cr a s). The first three columns are reproduced at the beginning of the three parts of this table. Each 
sub-table concerns one of the three EPIC instruments and follows the same arrangement. Cols. 4-16 (resp. 17-29 and 30-42) give the equivalent logarithmic likelihood L2 for 
the given instrument, the total count rate cr in the whole energy band (0.5-10.0keV) and its associated error (cr), the count rates in the different energy bands (Sx: [0.5-1.0keV], 
Mx: [1.0-2. 5keV], Hx: [2.5-10.0keV]) and their errors and, finally, the two hardness ratios HR\ and HR2 (Eqs.^and|2ji as well as their related errors (cthri and <x WR2 )- The 
count rates and the related uncertainties are all expressed in 10~ 3 cnts -1 . Col. 43 provides the total equivalent logarithmic likelihood L^ IC that has been compared to the 
adopted detection limit (see Tabled according to the combination of instruments used (Col. 44). For those sources for which the fit is improved while adjusting an extended 
psf model, a note ('ext') is stated in Col. 45. Finally Col. 46 provides a warning for some sources for which no instrument combination were found to be far enough from gaps, 
instrument edges or detector bad columns. The electronic version of the table further provides two additional columns that give the derived source coordinates with a precision 
of S .01 in right ascension and of 071 in declination. The full table can be obtained from the Centre de Donnees astronomiques de Strasbourg (CDS, http://cdsweb.u-strasbg.fr). 
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Table 4. Same as Table[3]but for the extended sources. cr psf (Col. 3b) gives, in arcsec, the <x extent of the adjusted Gaussian. 
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note that this procedure does not modify the spatial response 
of the detectors and that the known variations of the EPIC in- 
strument sensitivity with the axial distance will of course still 
affect our results. 

(iii) In the few cases for which multi-source fitting was rele- 
vant, we adopted the results obtained with this fitting. We how- 
ever paid a special attention to reject cases of fake multi-fitting 
sometimes induced by near-gap/edge effects or by multiple en- 
tries for a unique X-ray source in the preliminary source list. 

(iv) We finally checked every source in the final list by indi- 
vidually looking at the different image combinations. We elim- 
inated the very few double entries in the list. Doing this, we 
noticed a couple of presumably physical sources that were ig- 
nored by the detection algorithm. We decided to include those 
sources in the input source list of the emldetect task. Most of 
them were satisfactorily fitted, giving an equivalent logarithmic 
likelihood above the adopted detection threshold. These addi- 
tional sources were included in the final catalogue. 

(v) The main X-ray catalogue presented in Tablets based 
on the point-like source detection only. For some sources, the 
equivalent logarithmic likelihood L 2 is significantly increased 
if one adjusts an extended source model rather than a point-like 
model. These sources are flagged in Tableland we provide, 
in Table |4j a complementary extended source catalogue that 
gives, in addition to the results listed in the main catalogue, the 
emldetect extended-psf fit results for these sources. 

The final catalogue lists 610 sources in the XMM-Newton 
FOV, among which 19 are flagged as extended. Based on the 
edetectjzhain results, it provides, among other information, the 
source position, the total count rates in the different instruments 
and the two hardness ratios : 



HR, 



HR 2 = 



M- 



x 



> x 



M x + S x 
H x -M x 



H x +M 



(1) 
(2) 



x 



A sample of the catalogue is provided in Table [3] while 
Table 0] gives the complementary catalogue for the 19 ex- 
tended X-ray sources detected. In addition, source X#234 
appears clearly double in the EPIC image though it is not 
detected as an extended object. Table [5] is available on- 
line via the Centre de Donnees astronomiques de Strasbourg 
(CDS, http://cdsweb.u-strasbg.fr). Finding charts for the X-ray 
sources are provided by Figs.|2]and[3] 

4.2. Source Identification 

To determine the optical counterparts of the detected X- 
ray sources, we cross-correlated our source list with several 
existing optical/infrared c atalog ues. We used the US Naval 
Observatory dMonet et alJ 120031 USNO B1.0), the 2MASS 
All Sky Data Release dCutri et al.ll2003l) and the Guide Star 
Catalogue-II JGSC2.2||2001|). We also made use of the optical 
catalogue of Sung et alJ i 19981 SBL98 hereafter). However, 
the star positions in the SBL98 catalogue as available from the 
Centre de Donnees astronomiques de Strasbourg (CDS) show 




1 2 3 4 5 6 7 
r (") 
USNO B1.0 



10 01234567B9 10 

r (") 
SSB06: V<30 




Fig. 4. Cumulative distributions (solid lines) of the number (AO 
of closest associated counterparts as a function of the separa- 
tion radius (r) and for the different catalogues used. The hor- 
izontal dotted lines show the number N x of X-ray sources in 
the catalogue. The dashed lines, from top to bottom in each 
panel, correspond respectively to the best-fit <t>(r) function, to 
the number of truly associated counterparts <t> tlue and to the 
number of spurious ones <J> spur ., as a function of the correla- 
tion radius r. Dotted vertical lines show the correlation radius 
adopted for the purpose of source identification. 



clear shifts compared to the true positions on the sky. This re- 
sults from an excessive rounding of the star coordinates in the 
CDS database: they are given with a precision of respectively 
one second and one tenth of arcmin on the right ascension and 
declination. This is far insufficient in such a crowd ed field as 
NGC 6231. We therefore used the original SBL98 catalogue, 
obtained from the authors and that lists object coordinates a 
hundred times more pre cisely. Beyond the 860 objects with 
V < 16 listed in lSBL98l this UBV(RI) C & Her catalogue was 
comple ted with 7199 objects, extending the first version o f the 
ISBL98I catalogue down to V = 21. However, the ISBL98I field 
of view was limited to a 20' x 20' area and thus does not cover 
the whole EPIC FOV. It can thus not be used for identification 
throughout the entire field and we selected the X-ray sources 
that are located in the sub-region of the FOV that is covered by 
the SBL98 v2 catalogue. This yields a number of X-ray sources 
A^ BL = 536 as quoted in Table|5] More recently, one of us (H. 
Sung) acquired new UB V(RI)c observations covering a field of 
about 40' x 40' around NGC 6231. 30866 stars were observed 
down to V < 22. These observations will be presented in a 
forthcoming paper (Sung et al. 2006, in preparation) and we 
only focused here on the resulting catalogue. We will refer to 
this new catalogue as the SSB06. 

As a first approach, we investigated the possibility of sys- 
tematic differences between the reference frames of the differ- 
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Fig. 2. Detected X-ray sources overlaid on a DSS 1 image (V band) of the XMM-Newton FOV. The sources are indicated by black 
circles, with a radius of 3", similar to the adopted cross-correlation radius. The numbers above these circles give the internal 
X-ray source identification as provided in the first column of Table[3] The three blue circles indicate regions with a radius of 5', 
10' and 15' around X#279 (HD 152248). North is up and East is to the left. 



ent catalogues. For this purpose, we selected the bright O-type 
stars in the different source lists and we compared their loca- 
tions to the ones of their X-ray counterparts. Neither a signifi- 
cant systematic shift nor a field rotation was apparent. Typical 
1-cr dispersions computed on the differences between the lo- 
cations of the X-ray sources and their optical counterparts are 
about 0.9" in right ascension and 0.7" in declination. Similarly, 
the 1-cr dispersion on the field rotation is about 3'. As a second 
step and for each of the previously cited catalogues, we deter- 
mined the closest optical counterpart of each X-ray source in 



the field of view. We then calculated the cumulative distribution 
(O(r)) of the closest associated counterparts as a func tion of the 
individual correlation radius (see ljeffries et alJl997h. The gen- 
erated diagrams are shown in Fig. HI Following Ijeffries et all 
we assumed that €>(>) is formed by two terms: the cumulative 
distribution of true correlations <& tme and the cumulative num- 
ber of spurious associations <t> S p U1 -.. This is expressed in the sim- 
ple relation: 

$0) = O trae + <D spul , (3) 
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Fig. 3. Same as Fig.|2] zoomed on the inner part of the FOV. The circle radius is 5'. From left to right, the three crosses respectively 
give the position of the geometrical centre of the cluster, its X-ray emission centre (computed adopting the pn-equivalent count 
rates for each source), and the position of HD 152248 (X# 279). 



it fits reasonably well the rising branch of the different curves 
plotted in Fig.0] Table|5]gives the obtained values of the A, B 
and <x parameters. As the hypothesis of constant B throughout 
the FOV is clearly violated in the case of NGC 6231, we also 
estimated the number of spurious associations using a more 
empirical approach. We arbitrarily shifted the X-ray source po- 
sitions by 30" in any given direction and we re-ran the cross- 
correlation at a fixed r con (either 2'.'5 or 3'.'0 according to the 
catalogue considered). The obtained number of spurious asso- 



J +(jV x -A)[l-exp(-7rr 2 B)] (4) 



l-exp|^ 



that can be adjusted to the empirical distribution. In Eq.@] Ax 
is the number of X-ray sources while A is the number of true 
correlations with the optical/infrared catalogue. B is the opti- 
cal/infrared catalogue density and cr is related to the statistical 
uncertainty on the X-ray source position. Though Eq.0]is ap- 
proximative and rests on the hypothesis of a uniform optical 
population (i.e. constant B and constant psf over the full FOV), 
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Table 5. Best fit parameters (Cols. 3-5) of the <5(r) function (Eq.0Ji for different optical/infrared catalogues. The catalogue name 
is given in Col. 1 along with the relevant number of X-ray sources in the corresponding field (Col. 2). The adopted correlation 
radius (r con ) for identification is given in Col. 6. Col. 7 lists the actual number N con of identified X-ray sources (see Table [6} 
and the corresponding percentage related to the considered number of X-ray sources (A^x)- Col. 8 gives the number of associated 
counterparts predicted by the distribution O(r) (Eq.|4} at a radius equal to r corr . It also provides the corresponding percentage of 
theoretically identified X-ray sources. The next two columns provide, among the number of associated counterparts 3>(r colr ), the 
predicted number of true (tf'trueC'corr)) and spurious ((DspurXrcorr)) counterparts. The contribution of true and spurious counterparts 
to the total number of (theoretically) associated optical sources are also given in the corresponding columns. 
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2MASS 
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0.91 
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384 
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383.8 (62.9%) 
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90.5 (23.6%) 
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384.3 
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3.5 
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383.7 (62.9%) 


362.7 (94.5%) 


21.1 (5.5%) 
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383.1 
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0.8 
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344 
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367.8 (60.3%) 


362.8 (98.6%) 


5.0(1.4%) 


SBL98 v2 


536 


431.8 


1.09 


11.9 


3 


447 


(83.4%) 


451.9 (84.3%) 


422.1 (93.4%) 


29.8 (6.6%) 


SSB06: V < 19 


610 


396.6 


1.01 


4.3 


2.5 


384 


(63.0%) 


395.7 (64.9%) 


378.4 (95.6%) 


17.3 (4.4%) 


SSB06: V < 20 


610 


422.6 


0.95 


13.9 


2.5 


450 


(73.8%) 


453.7 (74.4%) 


408.9 (90.1%) 


44.9 (9.9%) 



ciations is n ever larger by more than 10% than the one esti- 
mated by the Jeffrie s et alJ (119971) method . 

The SSB06 catalogue is too dense for the relative crudeness 
of the X-ray source positions (cr a6 = Of.' 7 + Of.' 3 on average, cr aS 
being defined on 2D). Indeed, even adopting a limited cross- 
correlation radius of 2'.' 5 would yield over 100 spurious identi- 
fications. We thus decided to decrease the limiting magnitude 
of the catalogue. The maximum of the true function is obtained 
adopting V < 20. At the distance of the cluster, this corre- 
sponds to the magnitude of a M0 dwarf (M ~ 0.5 M G ). PMS 
low-mass stars being brighter than ZAMS stars of the same 
mass, the progenitors of M0 stars should thus still be included 
in the optical list. Finally, we note that a significant improve- 
ment (in terms of the relative percentage of spurious associa- 
tions) is obtained when restricting the SSB06 cross-correlation 
to objects with V < 19. The drawback is that the number of 
true associations is also significantly reduced. Table[3]lists the 
best fit parameters of Eq.|4]forboth cases and Table[3]provides 
the SSB06 cross-identifications down to V < 20. We leave to 
the user the choice to restrict the list to V < 19 according to 
his/her motivations. 

From the cumulative distributions shown in Fig. @J we 
adopted the cross-correlation radii corresponding to the knees 
in the distributions of counterparts; these are reported in Table 
[5] The percentage of identified sources ranges from 55 to 83% 
according to the catalogue used. The results for the SBL98 v2 
and the SSB06 catalogues are clearly in contrast with the other 
catalogues. With about 75% of the total number of X-ray 
sources in the FOV being identified, among which less than 
10% statistically correspond to spurious associations, the lat- 
ter catalogue is probably the most appropriate for the identi- 
fication processes. In the following, we thus adopt the SSB06 
catalogue, that covers the complete EPIC FOV, as the main ref- 
erence in the identification of the sources. Table[6]provides the 
cross identifications between the X-ray source lists and the dif- 
ferent optical catalogues. We find at least one counterpart in 
one catalogue for about 85% of the X-ray sources. 

While carrying out this work, we noticed some confusion 
between the names of several sources reported in the widely 



consulted SIMBAD database. For this reason, Table [6] also 
gives other commonly adopted source denominations such as 
HD, CPD and Braes numbers. The Seggewiss numbering is 
also extensively used in the literature related to NGC 623 1 . We 
therefore used the original chart of Seggewiss ( 1968b) - subse- 
quently completed by Rabou det al J i 19971) - and we rederived 
the cross-correlation to avoid any previous misidentification. 

4.3. The detection limit 

This paragraph aims at the evaluation of the detection limit of 
the present X-ray catalogue. Though essential, this question is 
far from trivial because the detection limit is, a priori, not uni- 
form throughout the field of view. Besides the areas where the 
detectors do not overlap and the presence of gaps between the 
detector CCDs, the XMM-Newton effective exposure duration 
is decreasing from the FOV centre towards its edges. In addi- 
tion, both the crowdedness of the field in its central part and 
the numerous bright sources preferentially located in the core 
of the cluster also affect the detection limit in a non uniform 
way. As an approximation we neglect the effects of the gaps, 
mainly focusing on the three other effects. 

The exposure maps computed for the three EPIC instru- 
ments and their different combinations display a smooth de- 
crease from the centre of the detector to its edges by about 
a factor of three. In terms of the amount of signal collected 
for two similar sources - one located near the FOV axis, the 
other near one of its edges - the number of counts n received 
will be three times higher near the axis. Neglecting any back- 
ground effect, the signal-to-noise ratio is approximately given 
by S IN - V"- For the outer source, it is therefore smaller by 
a factor of V3 ~ 1.7. To the first order, the detection limit in 
the outer parts of the field is thus about a factor two higher 
than in the central part of the FOV. As a next step, we used the 
SAS task esensmap to build sensitivity maps corresponding to 
the current exposure maps and to the adopted logarithmic like- 
lihood detection thresholds L^. The sensitivity maps obtained 
actually provide the minimum number of counts for a source 
to be detected by the detection task emldetect according to the 



Table 6. Cross-identification of the X-ray source catalogue of Table |3] with the different optical/infrared catalogues: 2MASS, GSC2.2, USNO B1.0 and SSB06 (V < 20). d cc 
is the distance (in ") between the positions of the X-ray source and of its closest counterpart. Only the closest identification has been reported in this table. Cross-identification 
with commonly used denominations is given in the last columns (Col. 26 to 31) of the table. 
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Fig. 5. Bottom part of the distribution of the pn- equivalent 
count rates of the X-ray sources as a function of their distance 
(d) to HD 152248. The plain line shows the adopted lower limit 
given by Eq.|5] 



given equivalent logarithmic likelihood threshold. These maps 
indeed predict that the sensitivity of the EPIC camera is twice 
larger near the axis than in the outer parts of the detector what- 
ever the instrument combination is. This is in agreement with 
our previous estimate. 

Accounting for the variation of the source density and the 
distribution of the bright sources in the FOV is a more tricky 
issue. We chose to adopt a completely empirical approach, tak- 
ing advantage of the large number of X-ray sources in the field. 
We assumed that a very good indication of the detection limit in 
the different parts of the field is given by the brightness of the 
faintest sources detected in these selected areas. We adopted 
the following approach. Because of the presence of gaps, we 
computed an equivalent EPIC pn count rate, in the range 0.5- 
10.0 keV, for each source. To the first order, the relation be- 
tween the count rates measured in any of the two MOS de- 
tectors and in the pn detector is approximately linear. Using 
the count rates obtained for sources that were detected on sev- 
eral EPIC instruments, we thus derived an empirical conversion 
factor between the MOS1, MOS2 and pn count rates. These 
factors were then used to obtain the so-called pn-equivalent 
count rates for all sources and, in particular, for those that fall 
in the gaps of one or several instruments. This yields approxi- 
mately coherent count rates for the different sources, whatever 
their position on the detectors. Figure|5]displays the source pn- 
equivalent count rates as a function of the distance from the 
central axis of the FOV - assumed to be the position of the 
binary HD 152248. A lower limit is clearly seen in the distri- 
bution. Selecting the faintest sources (i.e. those displaying the 
lowest equivalent count rates) in successive rings centered on 
HD 152248 provides an approximate sampling of this limit. We 
then adjusted a polynomial and derived an empirical detection 



Fig. 6. Lower panel: Estimated detection limit expressed in 
terms of the observed flux (in log(erg cirT 2 s )). Upper panel: 
Equivalent detection limit, expressed in ISM-absorption cor- 
rected luminosity (in log(erg s -1 )), for sources located in the 
NGC 6231 cluster {DM = 11.07, « h ,ism = 0.26 x 10 22 cnr 2 ). 
The different lines refer to the different mekal model tempera- 
tures adopted for the conversion. The energy band considered 
in both panels is 0.5-10.0keV. 

limit in terms of pn-equivalent count rates (crr lm ) as a function 
of the distance (d) from the field axis. This limit (in units of 
10~ 3 cnt s _1 ) is described by the following relation: 

cn im .(d) = 2.49214 - 0.65577 c/ + 0.1 1822 d 2 

- 0.00929 d 3 + 0.00030 d 4 { ' 

where d is the distance to HD 152248 expressed in arcmin. 
Eq.|5]is shown in Fig. [5] Clearly the detection limit is higher in 
the central part of the field (d < 5'), most probably because of 
the higher source density and because bright sources are pref- 
erentially located in the inner part of the FOV. The sensitivity 
also decreases towards the CCD edges, as indicated both by 
the exposure maps and the sensitivity maps. Finally we used 
single temperature optically thin thermal plasma models of the 
Raymond-Smith type to convert the pn-equivalent count rates 
given by Eq. |5] to fluxes and luminosities. For this purpose, 
we adopted the conversion computed by the WebPIMMS con- 
verter 3 , assuming a column density of 0.26 x 10 22 crrT 2 , typical 
of the interstellar absorbing column for the cluster. Results are 
displayed in Fig. |6]for three different plasma temperatures. In 
conclusion, the flux detection limit is approximately located 
between 3 x 10~ 15 and 1.5 x 10~ 14 erg cm 2 s _1 , depending on 
the distance from the detector axis and on the source spectrum. 
In the central part of the FOV, we consider that the typical lim- 
iting flux is about 6 x 10~ 15 ergcrrT 2 s _1 for soft sources. 

3 WebPIMMS is a NASA's HEASARC tool powered by PIMMS 
v3.6a. It is hosted at the following URL: http://heasarc.gsfc.nasa.gov/ 
Tools/w3pimms.html 
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Fig. 7. 2-D distribution of the surface density of the X-ray 
sources. The image is centered on the location of HD 152248. 
The original source density chart was convolved with a 
Gaussian with <x = 1'. Overplotted contour levels correspond 
to 1, 1.5, 2, 2.5, 3 and 3.5 sources per arcmin 2 . 
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Fig. 8. Plain line: cumulative distribution of the number of X- 
ray sources (N) with increasing distance (d) from HD 152248. 
Dashed line: idem, but computed assuming a uniform distribu- 
tion of the 602 sources (dotted horizontal line) located within a 
15' (dotted vertical line) circle around HD 152248. 



5. The X-ray sources in NGC 6231 

This section presents an overview of the main properties of the 
X-ray sources in NGC 623 1 . No attempt will be made here to 
investigate the characteristics of the different sub-populations 
of the cluster. This work is postponed to future devoted papers. 



5.1. Spatial distribution 

As seen from Figs. ^ to [3] there is an obvious clustering of 
the X-ray sources in the inner part of the FOV. Their spatial 
distribution projected on the sky presents, at first sight, a rev- 
olution symmetry around the centre of the field, i.e. the po- 
sition of HD 152248. Considering the sources located at less 
than 15' from HD 152248, we computed the geometrical cen- 
tre of the source distribution. We also computed the brightness 
centre of the X-ray image. For this purpose, we adopted the 
pn- equivalent count rates for each source. The two centres are 
located slightly East of HD 152248, at no more than 30" (see 
Fig. |3} ■ From the two-dimensional map of the X-ray source 
density (Fig. 0, we conclude that there is only a slight de- 
viation from this scheme and that the X-ray source distribu- 
tion shows a slight N-S elongation. In the following, we how- 
ever assume that the distance from the cluster centre, i.e. from 
HD 152248, remains the main parameter that shapes the source 
distribution. We also adopt the position of HD 152248 as the 
very centre of the cluster. 

From Figs. [7] and [8j it is clear that the radial distribution of 
the sources is not uniform and that most of them lie within a 10' 
radius around the cluster centre. We computed the radial den- 
sity profile of the X-ray emitters and we adjusted an empirical 



King density profile dKindl 19621) for a spherically distributed 
source population : 

12 



fid) = k 



i + (d t /d c y 



(6) 



where k is the central density, d c the core radius and d t the 
limiting radius. The King profile is very sensitive to k and d c , 
but less sensitive to d t which is indeed less meaningful for 
open clusters in the Galactic plane. The best fit parameters are 
k = 8.9 arcmin -2 , d c = 6'5 and d t = 2Q'.5. As indicated by 
Figs. |5] and |6] our detection limit depends on the location of 
the source on the detector. In a second step, we thus applied a 
relative correction to the X-ray density profile, accounting for 
the sensitivity difference as a function of the distance to the de- 
tector axis (crosses in Fig. [5}- The profile is now sharper and is 
described by: k = 7.6 arcmin -2 , d c — 3'1 and a! t = 1.5x 10 3 ar- 
cmin. In Fig. [9] we also present the density profile of the stars 
in SSB06. Investigating the source density distribution as a 
function of the limiting magnitude of the catalogue and of the 
distance to the detector axis, we further note that the SSB06 
catalogue is almost undoubtedly incomplete in the field centre 
above V — 18. This is easily explained by the number of bright 
sources (V ~ 5 - 10, see Fig. |3J in this region, that renders 
the detection of faint sources more difficult. For this reason, 
Fig. [9] is restricted to objects brighter than 17 in the V band. 
NGC 623 1 is further embedded in the Sco OB 1 association. 
As a consequence, the surface density does not drop to zero 
in the outer region of the field. We thus subtracted a threshold 
of 2 arcmin -2 prior to the adjustment. King best-fit values are 
this time k — 8.6 arcmin -2 , d c =2'.7 and d t = 1.4 x 10 3 arcmin. 
From Fig. |5] the correlation between the X-ray and op- 
tical surface density profile is obvious and yields similar 
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Fig. 9. Surface density profiles of the X-ray sources (open cir- 
cles) and of the optical sources (open squares) with V < 17. 
Crosses indicate the X-ray density profile corrected for the em- 
pirical EPIC sensitivity curve (see text). Best-fit King profiles 
are overplotted. 



core radii for NGC 623 1 . It further suggests that most of the 
detected X-ray emitters are physically belonging to NGC 623 1 . 



As discussed in e.g. ISung et all d2004l) . X-ray emission 
is probably one of the best membership criterion for young 
stars in open clusters. The present X-ray observations proba- 
bly provide the best census so far of PMS stars in NGC 623 1 ; 
though this census is probably still incomplete. However, the 
NGC 623 1 X-ray sample might be contaminated by foreground 
(field stars) and background (AGNs) objects. As a last check, 
we thus roughly estimated the probable number of foreground 
and background X-ray sources detected in the present cam- 
paign. Starting with the foreground objects, we proceeded as 
explained below. Accounting for the different typical X-ray lu- 
minosities for field stars of spectral type O to M and for our flux 
detection limit, we estimated the maximum distance at which 
a star can be located while still being detected. Using the so- 
derived distance, we computed the volume projected onto the 

XMM-Newton FOV. As a last step, we adopted t ypical star d en- 

\ ii — 1/ 



sities in the solar neighbourhood as quoted by Allen ( 1973| for 
the different spectral types. We finally end up with about 20 
foreground X-ray sources, most of which are expected to be 
G-type objects (12 stars) and F-type dwarves (4 or 5 stars). 
However, the previous approach does not account for proba- 
ble active stars or RS CVn in the FOV, which have lower spa- 
tial densi t ies but much higher luminosities. Using the work of 
Makarov (2003b), we found that about 21 galactic active stars 
could be detected in the EPIC FOV. This yields a total of ap- 
proximatively 41 contaminating galactic sources. As an addi- 
tional check, we also used the X-ray stella r log jV(> S) - logS 
curve at low galactic latitudes provided bv lMotch et alJ J2003). 




-2 

log (cr eq (cnt s -1 )) 



-1 



Fig. 10. Cumulative distribution of the number of sources as a 
function of increasing pn-equivalent count rate. The horizontal 
dotted line indicates the total number of sources in the FOV. 

Again we found that about 40 galactic X-ray sources are to be 
expected within our EPIC FOV. 

We then obtained a rough estimate of the number of extra- 
galactic background sources in our EPIC field. The Galactic 
coordinates of the cluster are l a = 343?46, b u = + 1° 19. 
Therefore, the total neutral hydrogen column density along 
this direction must be extremely large and should produce 
a substantial absorption of X-ray photons from extragalactic 
sources. Although they are in principle not suited for lines of 
sight at \bn\ < 5°, we used the D IRBE/IRAS extinction maps 
provided by ISchleeel et alJ dl998h to derive a first order esti- 
mate of the total column density. In this way, we find a total 
Galac tic E(B - V) of ab out 5.6 mag. Using the gas to dust ra- 
tio of iBohlin et all dl978l) . we thus estimate a neutral hydro- 
gen column density of ~ 3.2 x 10 22 cirT 2 . Assuming that extra- 
galactic background sources have a power-law spectrum with a 
photon index of 1 .4, and suffer a total interstellar absorption of 
3.2 x 10 22 crrT 2 , the mean detection limit 1.9xl0 _3 cnts _I with 
the pn camera translates into unabsorbed fluxes of 1.2 x 10~ 14 
and 3.5 x 10~ 13 ergcnr 2 s" 1 in the 0.5 - 2.0 keV and 2.0 - 
10 keV band re s pectiv ely. Using the log N - log S relation from 
iGiacconi et alJ d200 ll) . one expects thus about 13-16 extra- 
galactic objects among the detected sources. Thus, about 2% 
of the total number of sources could be background AGNs. 
It should be emphasized that these background AGNs are ex- 
pected to appear as rather hard (i.e. heavily absorbed) X-ray 
sources. 

In summary, both the geometrical and X-ray brightness 
centres of the set of detected sources correspond to the optical 
cluster centre. The radial profile of the X-ray source density is 
well correlated with the optical source radial profile. Both indi- 
cate a cluster core radius close to 3'. Finally, we expect that less 
than 10% of the presently detected sources correspond to fore- 
ground or background objects. We thus propose that the large 
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majority of the X-ray emitters revealed by the present XMM- 
Newton campaign are mostly belonging to NGC 623 1 . Some of 
them might alternatively belong to the Sco OB 1 association, in 
which NGC 6231 is embedded. 



5.2. Emission distribution 

While HD 152248, the brightest X-ray emitter in the FOV, dis- 
plays a pn-equivalent count rate larger than 0.36 cnts -1 , most 
of the other sources are much fainter with a count rate be- 
low 10~ 2 cnt s -1 (Fig. [Tot. It is clear from Fig. ^ that most of 
the brightest sources - associated with the O-type objects of 
the cluster - are relatively soft while the majority of the X- 
ray emitters have their maximum of emission in the medium 
band. Except for the brightest sources, there is no obvious cor- 
relation between the source intensity and the source hardness 
ratios. On average, the detected sources are moderately hard 
with HR X > and HR 2 < (Fig. [TTJ- The hardness ratios 
might however show a slight increase towards the edge of the 
detectors, probably due to the relative dominance of low-mass 
stars in the outer regions of the FOV. The histograms of the de- 
tected source count rates in the Sx, Mx and Hx bands (Fig.1121 
reveal clear peaks around 0.7, 1.0 and 0.2 x 10~ 3 cnts~' re- 
spectively. The count rate in the 0.5-10.0 keV band clusters 
at 2 x 10~ 3 cnts" 1 and the two hardness ratios around 0.2 
and -0.6 respectively. Accounting for the cluster typical ISM 
absorbing column «h,ism = 0.26 x 10 22 cm -2 , these values 
are roughly described by a mekal model with a temperature 
of kT = 1.0 - 2.0 keV. The corresponding observed flux is 
about 5 x 10~ 15 ergcirT 2 s~'. Adopting a distance modulus 
DM - 11.07, this yields a luminosity log(L x ) ~ 30.5 (ergs -1 ) 
for a typical X-ray emitter in the cluster. 

6. Summary 

We presented the first results of an XMM-Newton campaign 
on the young open cluster NGC 623 1 in the Sco OB 1 associ- 
ation. With an effective cumulated exposure time of 175 ks in 
the two EPIC MOS instruments and of 147.5 ks in the EPIC 
pn, the campaign was split into six successive observations ac- 
quired within 5 days. The combined image, built from the data 
collected by the three EPIC instruments during the six point- 
ings, reveals an extremely crowded field. Running the S AS task 
emldetect, we detect 610 X-ray sources among which 19 are 
reported as extended. The latter are probably constituted by 
non-resolved point-like sources rather than by physically ex- 
tended sources. We present an X-ray catalogue covering the 
XMM-A^evvfon FOV and we cross-correlate it with several op- 
tical/infrared catalogues. We find at least one optical and/or 
infrared counterpart for more than 85% of the X-ray sources 
within a limited cross-correlation radius of 3" at maximum. We 
estimate our detection flux limit to lie between about 3 x 10~ 15 
and 1.5 x 10" 14 ergcm" 2 s _1 depending on the position on the 
detectors and on the source spectrum. 

The surface density distribution of the X-ray sources peaks 
at the centre of the cluster, which we find to be located 
very near HD 152248, and presents a slight N-S elongation. 
Concerning the radial profile of the surface density distribution, 
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Fig. 12. Distribution of the X-ray emitter count rates in the dif- 
ferent energy bands considered. The four panels were plotted 
only using the 392 sources detected with the pn instrument. 
Results for the EPIC MOS instruments are similar. The last bin 
includes the contributions of all the brightest X-ray sources. 



over 50% of the sources are confined within a 6' radius from 
the cluster centre and about 80% within 10'. The estimated con- 
tamination by foreground and background objects is about 9%. 
There is a good similarity between this radial profile and the 
distribution of stars brighter than V - 11, suggesting that most 
of the sources physically belong to NGC 623 1 . The radial sur- 
face density profile of the X-ray sources is well described by a 
King profile with a core radius of about 3', similar to the one 
indicated by the V < 17 optical source density profile. 

Finally, beside a few bright and soft objects correlated with 
the O-type stars of the cluster, the large majority of the X-ray 
population is relatively faint (pn-equivalent count rate below 
10~ 2 cnt s _1 ) and displays an intermediate spectrum of a typical 
temperature probably around 1.0-2.0keV. Typical count rates 
for the sources are 2.0, 0.7, 1.0 and 0.2 x 10~ 3 cnts -1 respec- 
tively in the total energy band (0.5-10.0 keV), and in the three 
energy sub-ranges 0.5-1.0, 1.0-2.5 and 2.5-10.0 keV. At the 
NGC 623 1 cluster distance, these values roughly correspond to 
an X-ray luminosity of about log(Lx) ~ 30.5 (erg s _1 ). 

More detailed investigations of the X-ray properties of the 
different sub-populations (early-type stars, PMS objects) of the 
cluster will be presented in subsequent papers in this series. 
Finally, the X-ray data related to specific early-type binary sys- 



tems of particular interest are (will be) presented in dedicated 



paper s (see e.g. lSana et al!2004t 1*2005 



)P 

a,2 



2006a as well as San 



2005), together with the derivation of their orbital and physi- 
cal parameters obtained on the basis of an extensive spectral 
monitoring campaign in the optical domain. 

Acknowledgements. It is a pleasure to thank Dr. W. Seggewiss for 
providing a copy of his original paper on the cluster photometry 
<Seggewissi ri968rA our system manager, A. Detal, for essential as- 



18 



H. Sana et al. : An XMM-Newton view of NGC 623 1 . I. 




-3 -2 -1 5 10 15 5 10 15 -1 -0.5 0.5 

Log <CR (ont s-')) d (') d (') HR, 



Fig. 11. Left panel: Hardness ratios versus count rate. Middle panel: Hardness ratios versus distance (d) from the FOV centre. 
Right panel: HR2 versus HR\. The three panels are built for the sources detected with the pn instrument. Results for the EPIC 
MOS instruments are similar. 
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Appendix A: On correcting the L 2 values in 
SASv 5.4.1 

The equivalent (or transformed) logarithmic likelihood L2 as- 
sociated with each source candidate detected by the SAS task 
emldetect (column DET_ML in the output file) is given by: 

L 2 = -ln(l-p(^,L')) (A.l) 

with 

i=n 

L'=Yj lh (A ' 2) 

where P is the incomplete Gamma function, v is the num- 
ber of degrees of freedom (d.o.f.) of the fit, n is the number of 
input images (i.e. the number of energy bands times the num- 
ber of instruments considered), and /,■ = C, /2 wit h C, be i ng the 
Cash statistic for image i, specially designed by Cash (Il979h 
for photon counting experiments. More insight into the physi- 
cal meaning of Ea. lA.ll will be given in the next section. In this 
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section we focus on the implemented patch for correcting 
values. 

Indeed the logarithmic likelihood is known to be er- 
roneous in SAS version v 5.4.1 and earlier versions (XMM- 
Newton News #29 - 1 l-Mar-2003). According to SAS 
Observation Report SASv5.4/8665 4 , the factor 2 in equation 
U - Ci/2 has been forgotten, leading to erroneous U and 
hence Li. Knowing the number of degrees of freedom v, it is 
a simple exercise to invert Eg. IA. H and to obtain values for L' . 
From Eq. IA.2I it is obvious that the corrected value for L' is 
L' con = L' /2, to be used in Eg. lA.ll to recover the corrected 
L™" value that can then be used for scientific analyses. 

For large values of L2 (L2 ^ 10 000) the numerical limits of 
classical compilers are however exceeded. Fortunately Eg. IA.ll 
tends to a linear relation between L' and L2 for large values 
and for a given v. The correction is therefore straightforward 
with L C ° 1T = L2/2. Though this bug was present at the time we 
analysed the data, this issue has been fixed later in SAS version 
v6.0. We checked our corrected L^™ values against SAS v6.0 
and found them in close agreement. 

Appendix B: On the choice of coherent detection 
thresholds using the transformed logarithmic 
likelihood L 2 

As it can be deduced from the previous section (App. IaV 
the logarithmic likelihood L2 is related to the probability that 
a detected source candidate could be explained by pure ran- 
dom Poissonian fluctuations (and zero count in the source). 
Computed for each source of the input list, it uses a combi- 
nation of the Cash statistic C, obtained for the different input 
images i. The Cash statistic C, actually results from a likelihood 
ratio test and obeys ax 2 distribution dCashl 19791) with 3 or 4 de- 
grees of freedom (i.e. the intensity, the X- and Y-coordinates of 
the source and, eventually, the extent of the source if allowed). 
Therefore any linear combination of n Ci, and hence any com- 
puted 2L' , also follows a^ 2 statistic with n + 2 or n + 3 d.o.f. In 
this sense, the transformed logarithmic likelihood L2 is indeed 
linked, through the simple relationship 

£ 2 = -ln(0, (B.l) 
where 

e = e(H = 1 "*&4 (B - 2) 

to the probability Q for a random Poissonian fluctuation to have 
caused such a high value of 2V = Yi]Z" Ci as the one observed. 
The equivalent logarithmic likelihood L2 will therefore be large 
if the observed source is likely not produced by a statistical 
fluctuation, and small otherwise. 

As a consequence, a threshold in L2 can in principle be 
adopted as a detection limit. However, as we show below, while 
the expression given in Eg. IB. 21 indeed takes into account the 
number v of d.o.f. to compute the Q probability and the sub- 
sequent value of L2, it does not allow a direct comparison be- 
tween L2 obtained with different numbers of input images. This 
statement is illustrated in the following due consideration. 

4 http://xmm.vilspa.esa.es/xmmhelp/ 



Table B.l. Illustration of consistently determined L2 thresholds 
(Col. 2) for the different instruments and instrument combina- 
tions reported in Col. 1 . The number of input images («) and 
corresponding degrees of freedom (v) are given in Cols. 3 and 
4. L (Col. 5) is linked to L 2 through Eq. IA.1I A given L' is 
also linked to other L' of this table through Eg. lA.2l (see text). 
We emphasize that adopting any of the L2 or L' presented in 
this table automatically determines the other values of L' and 
L2 reported here below. 



Instr. Comb. 


L 2 


n 


V 


V 


MOSl 


10.00 


3 


5 


13.75 


MOS2 


10.00 


3 


5 


13.75 


pn 


22.77 


3 


5 


27.51 


MOS1+MOS2 


19.25 


6 


8 


27.51 


MOSl+pn 


31.70 


6 


8 


41.26 


MOS2+pn 


31.70 


6 


8 


41.26 


MOSl+MOS2+pn 


40.86 


9 


11 


55.02 



Let us assume that we are dealing, for example, with 3 en- 
ergy bands and let us only consider point-like source fitting 
(parameter withextendedsource-'no'). For the purpose of the 
demonstration, let us adopt a uniform detection threshold, for 
any instrument or instrument combination, of L2 - 10. 

As a first step, let us deal with the source detection on the 
EPIC MOSl images. In this particular configuration, there are 
three input images (n = 3) that correspond to the three energy 
bands. From the inversion of Eg. IA. II with Zj^ 0S1 = 10 and 
v = 5, we obtain £mosi = 13.75, where L' MQSl is the sum of the 
/ MOS1 for each of the three input images as given by Eg. IA.2l 
i.e. 

!=3 

rl _ V ;MOSl 

L MOSl _ H 
i=l 

Now assuming that the two instruments MOS 1 and MOS2 are 
exactly identical, a detection threshold L^ OS2 = 10 similarly 
corresponds to L' MQS2 = F 0S2 = 13.75. 

In a next step, let us work with a combination of the two 
EPIC MOS instruments. Eguation lA.2l allows us to easily build 

the combined ^mosi+mos2 as ^ e sum °^ t ^ le ^ ^ or eacn mstm ~ 
ment and energy band : 

j 1 _ \f /MOSl , V 1 /MOS2 _ 97 <ri 

MV10S1+MOS2 ~ / /i Z-i ' 

i=l i=l 

With two instruments and hence 6 images, ^mosi+mos2 
lows a x 2 distribution with 8 d.o.f. (v = 8). Equation IA.1I 
then gives /MOS1+MOS2 _ 19 25 q U j te different from the value 
£MOSi+MOS2 _ 10 obtained with the adopted constant threshold 
limit L 2 = 10. 

If we consider the use of two identical detectors, the fact, 
on one side, to combine them and, on the other side, to adopt 
the same statistical limit for both an isolated detector and a 
pair of them, allows us to go deeper. Actually, the combined 
logarithmic likelihood is twice the individual ones: 

^MOSl+MOS2 = 2^MOSl = ^MOS2- 

Thus, this kind of threshold does not preserve the detection 
limit which is dependent on the particular combination used. 



H. Sana et al. : An XMM-Newton view of NGC 623 1 . I. 



21 



If we want to preserve the detection limit adopted for a 
single instrument, we must, in this example, also multiply 
the detection threshold by a factor of two, adopting the value 
27.51 instead of 13.75 and consequently 19.25 instead of 10 
for the transformed L2 statistic. We can of course extend this 
result to the pn detector. Making the reasonable assumption 
that L' pn * 2L' MQS , a similar reasoning gives ZJJ" - 22.77, 

L MOSl +P „ = L MOS2 + p„ = 7Q L MOSl + MOS2 + p„ = ^ ^ 

far from the value of 10.0 initially adopted. The intermediate 
results and numbers of d.o.f. used in establishing these values 
are given in Table IB. II Basically, when combining several 
instruments together, we improve the Poissonian statistics. 
The fact of adopting a constant value for L2 for various 
instrumental combinations implies a cut-off in fluxes or count 
rates that is dependent on the number of detectors considered. 
Instead, if we prefer to stabilize the cut-off in absolute values 
of the signal rather independently of the combination used, we 
have to adapt the L2 value to the situation. 

In summary, one of the main results of the present discus- 
sion is that one should not adopt a constant threshold limit in 
L2 for different instrument combinations if one wishes to pre- 
serve the uniformity of a cut-off level adopted for a given in- 
strument or combination. We have shown that the L2 thresholds 
in different combinations are linked through Eas. IB.ll and lB.2l 
and through the detector physical characteristics that condition 
the C, values. In consequence, adopting a particular value as a 
threshold for a specific instrument or instrument combination 
implicitly assigns related values to the L2 thresholds for any 
other instrument or combination considered. Therefore, if one 
wants to adopt a consistent detection threshold whatever the 
considered instrument or combination are, the previous reason- 
ing becomes a forced step. This issue is particularly relevant to 
consistently deal with sources that fall on gaps or on specific 
detector areas where the different instruments do not overlap. 
We remind that this does not modify the spatial response of the 
detectors (nor the effect of the field crowdedness). Thus, spa- 
tial variations in the effective count rate threshold are still to be 
expected and, indeed, they are observed (see Figs. |5]and|6}. 

We finally remind the reader that the above presented 
method to determine self-consistent L2 thresholds rests on two 
simplifying, but reasonable, assumptions. The first is that the 
two EPIC MOS instruments are identical. The second is that 
the EPIC pn yields approximately L' pn * 2LJ^ 0S . Any refine- 
ment of these two assumptions (i.e. any relation giving the L' 
of one instrument as a more realistic function of the L' of the 
other instruments) can be easily included in the method. This is 
however beyond the scope of the present discussion. The pro- 
cedure illustrated here has been used in the making of our cata- 
logue. The figures appearing in Table|2]were indeed established 
in a similar way (adopting L^ osl = l¥° S2 = 11) and represent 
the threshold actually utilized. 



